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Abstract 

The last few years have seen a revolution in very -high y-ray astronomy (VHE; 
E> 100 GeV) driven largely by a new generation of Cherenkov telescopes (namely 
the H.E.S.S. telescope array, the MAGIC and MAGIC-II large telescopes and the 
VERITAS telescope array). The Cherenkov Telescope Array (CTA) project fore- 
sees a factor of 5 to 10 improvement in sensitivity above 0.1 TeV, extending the 
accessible energy range to higher energies up to 100 TeV, in the Galactic cut-off 
regime, and down to a few tens GeV, covering the VHE photon spectrum with 
good energy and angular resolution. As a result of the fast development of the 
VHE field, the number of pulsar wind nebulae (PWNe) detected has increased 
from one PWN in the early '90s to more than two dozen firm candidates today. 
Also, the low energy threshold achieved and good sensitivity at TeV energies has 
resulted in the detection of pulsed emission from the Crab Pulsar (or its close en- 
vironment) opening new and exiting expectations about the pulsed spectra of the 
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high energy pulsars powering PWNe. Here we discuss the physics goals we aim 
to achieve with CTA on pulsar and PWNe physics evaluating the response of the 
instrument for different configurations. 

Keywords: 



1. Pulsars and Pulsar Wind Nebulae at Very High Energy 

Pulsars and their synchrotron nebula have been extensively observed in ra- 
dio and X-ray wavelengths. The non-thermal synchrotron emission detected from 
these systems provides evidence for particle (lepton) acceleration to extremely 
high energies, susceptible to inverse-Compton (IC) scattering to target soft pho- 
tons in the surrounding medium and produce VHE y-ray emission. Recent ob- 
servations of Galactic sources with Imaging Cherenkov Atmospheric telescopes 
(lACTs) have revealed PWNe to be the most effective Galactic objects for the 
production of VHE y-ray emission. As recently as 2004, only the Crab PWN was 
detected with a steady y-ray flux above 1 TeV of (2.1±0.1stat)xl0^^cm^^s^^ 



11161 150n . The development of new sensitive lACTs in the last years has raised 
the number of likely PWNe detected to at least 27 sources, whereas many of the 
unidentified y-ray sources are widely believed to be PWNe (or old relic PWNe). 

These VHE PWNe are believed to be related to the young and energetic pul- 
sars that power highly magnetized nebulae (a few /iG to a few hundred /iG). 
In this scenario, particles are accelerated to VHE along their expansion into the 
pulsar surroundings, or at the shocks produced in collisions of the winds with 
the surrounding medium. Leptons accelerated to relativistic speeds interact with 
the low energy radiations fields produced by synchrotron, thermal or microwave- 
background origins. As a result, non-thermal radiation is produced from the low- 
est possible energies up to ~100 TeV. For a few /iG magnetic fields, freshly in- 
jected electrons (and positrons) create a small synchrotron nebula around the pul- 
sar which should be visible in X-rays, in contrast to an often much larger TeV 



nebula, generated by IC processes (for a recent review see [|34l1 ). Typically only 
young pulsars with large spin-down power (£'>10-^-^erg s^^) produce prominent 
PWNe fisfl. This is illustrated in Fig. [T] where the spin-down power of radio 
pulsars from the ATNF catalog [|40n is shown versus their characteristic age. Pul- 
sars powering VHE y-ray PWNe are marked with red bullets whereas pulsars in 
which periodic emission originating in the pulsar magnetosphere (as detected by 
the Fermi-LAT satellite, 100 MeV<E< 100 GeV) are shown in blue. 
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Figure 1: Spin-down energy loss rate of radio pulsars listed in the ATNF catalog versus their 
characteristic age. The blue dots overlaid indicate the pulsars detected at high energies with the 
Fermi-LAT telescope whereas the red ones mark the pulsars associated to TeV PWNe. 



The discovery of this large population of TeV PWNe has provided a large in- 
put for multi-wavelength models of particle evolution and acceleration in these 
objects. Theoretical MHD models (carried out by e.g. pllj) describe well the 
observed emission by a relativistic wind of particles driven by a central pulsar 
blowing into the ambient medium creating a termination shock. The particles are 
believed to be accelerated somewhere between the light-cylinder and the termi- 
nation shock, although this question is still open. The broadband spectrum of a 
PWN thus provides constraints on the integrated energy injected by the pulsar as 
well as on the effects of adiabatic expansion and the evolution of the magnetic 
field. The spectral energy distribution (SED) consists of two components, 1) syn- 
chrotron emission extending from the radio into the X-ray and, in some cases, the 
MeV band, and 2) IC emission producing GeV and TeV photons. The combined 
observations at VHE, X-rays and radio wavelength are crucial to constrain the 
evolution of the nebula magnetic field as well as the magnetic-to-kinetic energy 
conversion ll29ll . 

Pulsars are sources of non-thermal radiation which in some cases extends over 
many decades in energy, from radio bands to y-rays. Thanks to Fermi LAT obser- 
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vations pulsars are now acknowledged as a distinct class of galactic y-ray sources. 
Fermi LAT Second Source Catalog [T] contains nearly 90 individual pulsars in- 
cluding 27 firm detections of millisecond pulsars (MSPs). Different models of 
pulsar activity give different predictions for emission characteristics in terms of 
directionality and energy spectra. Due to the highly anisotropic (beamed) and 
energy-dependent structure of pulsar emission the observed radiation patterns, i.e. 
the light curves and spectra, depend strongly on the line-of-sight with respect 
to the spin axis. In most cases, a simple power-law with exponential cutoff lo- 
cated at a few GeV fits the LAT spectra satisfactorily. The y-ray light curves are 
mostly either single-peaked or double-peaked; their absolute phasing with radio 
light curves provide important clues to determine topology and spatial extent of 
the emission region. It came as a surprise that in most cases of MSPs their y-ray 
properties are similar to those of normal pulsars despite much weaker magnetic 
field strengths of the former. Moreover, unpulsed y-rays from 1 1 globular clusters 
have been detected; this is interpreted as due to the ensembles of MSPs in these 
clusters. Some support for such interpretation comes from the case of J1823- 
3021 A in globular cluster NGC6624, the first MSP in a GC detected so far as a 



source of pulsed y-rays [1330 • 

Little is known, however, about the properties of pulsars above ^10 GeV 
where LAT performance is limited. Recent results from VERITAS [23] and MAGIC 
[|20[ [21I1 have shown, for the Crab pulsar, a clear detection, up to 400 GeV of 
pulsed power-law type emission, as opposed to an exponential cutoff as might 
have been expected. 

1.1. PWNe physics goals with CTA 

Observations at VHE have revealed PWNe as the most efficient type of source 
in accelerating particles and producing y-ray through IC mechanisms, even al- 
lowing the detection of such systems outside our own Galaxy [39]. Moreover, 
the VHE part of the SED or IC peak contains information of the total energy re- 
leased by the pulsar in its life time. This is due to the slower cooling time of TeV 

—2 —1/2 

electrons (t ~(4.8 kyr)B_5E-pgY ) compared with the live-time of keV-emitting 

—3/2 —1/2 

electrons (t ^(\.2 kyr)B^' Ej^^y ). PWNe thus act like calorimeters at VHE 
and observations with Cherenkov telescopes, and in particular with CTA, permit 
the investigation of the time-evolution of such systems. This also implies a bet- 
ter understanding of the still unknown mechanism of how the pulsar releases its 
rotational energy from a cold relativistic wind blowing out the light cylinder to a 
kinetic one with Lorenz factor of Ri3xlO^ at the termination shock. Within the 



4 



context of temporal evolution of PWNe, two types of object seem to be emerging: 
IVyoung plerions such as the Crab nebula (3], SNRGO.9+0.1 H], MSH 15-52 
lll5n and the newly discovered Crab-like VHE y-ray sources SNRG2 1.5-0.9 and 
Kes 75 ll46ll and 2) evolved (extended and resolved) systems (i.e. with characteris- 
tic ages T>104 yr), such as Vela X (3, HESS J1825-137 U^, HESS J1718-385 
@] and HESS Jl 809-193 llioll. 

The young systems show in general a good match with the morphology seen 
in X-rays, while in the latter group, very often the pulsar powering the TeV ple- 
rion is found offset with respect to the center of the TeV emission, with large 
size ratios between the X-ray and VHE y-ray emission regions. The evolution of 
the supernova remnant (SNR) blast wave into an inhomogeneous ISM and/or the 
high velocity of the pulsar, together with a low magnetic field value ( ~5/iG), 
may explain these large offsets as being the relic nebulae from the past history of 
the pulsar wind inside its host SNR. Electrons responsible of the VHE radiation 

1/2 

(with energy Ee=(18 TeV)Erpgy to inverse Compton scatter CMBR seed photons 
to energies ExeV, to be compared with the required electron energy in a trans- 
verse magnetic field of strength B = lO^B 5 G, to radiate synchrotron photons 

— 1/2 1/2 

of mean energy E^eV^ Ee=(70 TeV)B g' Ej,gy) do not suffer from severe radiative 
losses and the majority of them may survive from (and hence probe) early epochs 
of the PWNe evolution. This interpretation has been further supported by the dis- 
covery of the spectral softening of the VHE nebula HESS J 1825-1 37 as a function 
of the distance from the pulsar nlW . All these new observations have allowed a 
fast development of the field and a better understanding of the physics of PWNe as 
well as opened new questions about e.g. particle transport and dynamics, energy 
loses, inter-play between the high spin-down energy pulsar and its surrounding 
nebula and interaction between PWNe and the hosting SNR. 

As the largest and most efficient population of VHE 7-rays sources, it is clear 
that a deeper understanding of PWNe is one of the key goals of the CTA observa- 
tory. The following research lines have been pursued and evaluated for the three 
different configurations described previously in this volume to asses the CTA best 
response: 



• Maximum and minimum resolvable size of the PWN. 

The majority of VHE PWNe have a very large size up to 1.2° depending 
on their evolutionary stage and proximity (see e.g. VelaX or HESS J1825- 
137). One of the key factors for the exponential increase of the PWN num- 
ber in the last years has been the large field of view of the third generation of 
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Cherenkov Telescopes (~5''), which has allowed the detection of the large 

y-ray halo produced by cooled electrons. Similar features will be required 
for CTA to be able to image the whole plerionic emission. At least a flat 
response up to 1 .5° should be aimed for. 

Together with a large field of view, to be able to resolve sub-structures on 
the VHE emission and thus sites of VHE y-ray production, a good angular 
resolution is required of at least <3'. Given the smaller size and complexity 

of the nebulae detected in X-rays, an angular resolution comparable to sys- 
tematic pointing error (5" per axis) should be achieved for CTA (as well for 
the study of composite SNRs, see below). 

• Electron cooling effect The large energy range covered by CTA will be cru- 
cial to understand cooling effects in the electron parent populations, includ- 
ing the resolution of internal structures and the possibility to disentangle 
between synchrotron or adiabatic losses. 

• Composite SNRs One of the most intriguing aspects of the new observa- 
tions is the interaction between the host SNR and its PWN, and how much 
of the 7-ray emission is originated by the SNR. A good angular resolution 
is required to understand those composed systems, which often are unre- 
solved when observed with the present instruments (i.e. Kes 75 with a shell 
radius of 1 .75' or SNR G21 .5-0.9, with a radius of 2.5', to compare with the 
current Cherenkov telescopes angular resolution of ~3'). 

• Population studies PWNe have been proven to be the most efficient Galac- 
tic sources in producing VHE y-ray. This large efficiencies and large size 
(detection sensitivity is proportional to 1/d, where d is the source size, for 
extended sources and not 1/d^) permit the detection of PWNe up to 50 kpc 
(i.e. PSR J0537-6910 on the LMC). Observations with CTA will allow an 
homogeneous sampling of the Galaxy. 

• PWN Modeling The new observations and results will constrain and pro- 
vide input for MHD simulations of PWNe, such as estimation of mag- 
netic field, electron population behind the 7-ray (and synchrotron) emission, 
spectral characteristics, etc . 

1.2. Pulsars physics goals with CTA 

One of the key issues in pulsar studies is the question on how far the non- 
thermal spectra of pulsed radiation from known y-ray pulsars actually extend in 
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the energy domain and what are the shapes of their cutoffs. The shape and posi- 
tion of the spectral cutoff vary depending on the emission model considered and 
the observer orientation. In polar cap models (e.g. [30]) electrons are acceler- 
ated at low altitudes and radiate y-rays via synchrotron-curvature radiation. Since 
these 7-rays are created in super-strong magnetic fields, magnetic pair produc- 
tion is unavoidable, and hence, only those secondary photons which survive pair 
creation (a few GeV for typical pulsar magneto sphere) escape to infinity as an ob- 
served pulsed emission. A natural consequence of the polar cap process is a super- 
exponential cutoff of the spectrum above a characteristic energy Eq. The way to 
maintain the magnetic pair creation above the gap and avoid super-exponential 
cutoffs is to extend the gaps spatially. This is a property of the slot-gap model by 



14111 . In outer gap models, e.g. [35], y-ray production is expected to occur fur- 
ther away from the stellar surface. In this case the cutoff is determined either by 
photon-photon pair production or by maximal energy reached via inverse Comp- 
ton scattering. Apparent absence of super-exponential cutoffs in the Fermi LAT 
pulsar spectra have called into question, therefore, the classical picture of polar 
cap activity in pulsars. 

One of the main motivations of CTA is to explore the energy range of 10 
GeV to 100 GeV, by using very large reflectors to decrease the energy threshold 
to a few tens of GeV. CTA should thus offer a unique opportunity to explore the 
most extreme energetic processes in pulsar magnetospheres and beyond their light 
cylinders. The goal is to formulate self-consistent electrodynamic theory of pulsar 
activity and the VHE domain should play an important role in it. 

1.2.1. Pulsed VHE from the Crab pulsar - the theoretical context 

Pulsed VHE emission from the Crab pulsar is of particular interest for pulsar 
physics. In the context of outer gap models its presence and properties depend 
strongly on the actual location of outer boundary of the outer gap with respect 
to the light cylinder. Above a few GeV it is synchro-curvature emission due to 



primary electrons which dominates all other intrinsic spectral components ||35, 



3611 ■ However, it is strongly attenuated by the dense field of magneto spheric soft 
photons. The emerging radiation is then dominated by IC scattering of secondary 
and tertiary pairs with the soft photon field. 

There have been other suggestions about possible sources of pulsed HE and/or 



VHE radiation from pulsars, like the slot-gap model ( 
with a force-free structure [|24ll . the striped wind model 



4111 ). the magnetospheres 
43I1 or IC scattering of soft 



thermal and non-thermal photons from the pulsar with ultra-relativistic electrons 
in its pulsar wind zone 



aipnotc 
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In all present-day outer-gap models the location of the outer boundary (as well 
as the inner boundary) of the outer gap is not determined self-consistently; there- 
fore, it is a free parameter. The particular choice of the outer boundary location 
results in a particular extent and shape of the HE- VHE radiation: The outer bound- 
ary located at around 75% of the light cylinder radiusjresults in the spectrum with 
exponential cutoff before even reaching 25 GeV i35ll . But the outer boundary 
closer to the light cylinder leads to a lower attenuation probability of the intrinsic 
VHE photons in the field of soft photons. In particular, the VHE spectral compo- 
nent in pulsed radiation from the Crab pulsar as discovered by VERITAS [23] and 
MAGIC HIH may be accommodated by shifting the outer boundary to a higher 



value (K. Hirotani, private communication); moreover, the cutoff shape changes 



from an exponential-like to a double-power-law (Fig. 15 in 112011 ). This does not 
mean, however, that a power law combined with an exponential cutoff to describe 
the pulsar spectra as inferred from Fermi LAT are invalidated with the VERITAS 
discovery. The power law with exponential cutoff is just one of several template 
shapes used in the spectral analysis of Fermi LAT pulsar data. A statement that the 
best fit for a given pulsar data is a power law with exponential cutoff does not rule 
out significant deviations from this shape, especially at the highest energy bins. 
With large error bars at the highest energy bins there may be plenty of room to 
accommodate an additional component (e. g. expected due to some ICS at VHE). 

The discovery of pulsed VHE emission shows that the outer gap models should 
undergo major modifications if the gaps are to be closer to the light cylinder than 
assumed before. More accurate treatment of the electrodynamics, including cur- 
rents and possibly going beyond the light cylinder is necessary. High quality spec- 
tral, phase-resolved properties of the VHE radiation from Crab will be essential 
to help to develop realistic models of the magneto spheric gaps. CTA is expected, 
therefore, to contribute significantly to the progress in our understanding of the 
magneto spheric activity of high- spin-down pulsars. 

1.2.2. Unpulsed VHE from Globular Clusters 

Globular Clusters are known to contain many MSPs and they are expected to 



be the sources of unpulsed VHE radiation e.g. 1251 14811. The most attractive glob- 
ular cluster harboring probably as many as 200 MSPs (more than 20 radio-pulsars 
are known to date) is GC 47 Tucanae (47 Tuc). The upper limit from H.E.S.S. on 
the VHE flux above 800 GeV from 47 Tuc is at the level 6.7 x lO-^m-^s"! fl. 
The unpulsed VHE flux expected from 47 Tuc in the model by [48] is the sum 
of two IC scattering components due to primary electrons up-scattering ambient 
starlight and CMB photons but its level depends heavily on the strength of ambient 
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magnetic field in the cluster. The H.E.S.S. UL appears stringent enough to con- 
strain the parameter space (number of MSPs vs. ambient magnetic field strength) 
in the simple model of 47 Tuc [48]. H.E.S.S. reported recently on the discovery of 
the VHE radiation from vicinity of Terzan 5 (Ter 5) [5] . This radiation may be due 
to MSPs in the cluster, however a significant offset from the location of Ter 5 has 
no satisfactory explanation so far. For other globular clusters, with lower content 
of MSPs, the expected flux level between 100 GeV and 10 TeV may not be within 
the reach of present-day lACTs. However, CTA may be sensitive enough to detect 
some globular clusters and even to resolve their spectral shapes at VHE. 
CTA will, therefore, in principle allow the detailed study of the major ingredients 
of globular clusters which link MSPs to the expected VHE emission. 



2. CTA response simulations 

2.1. Pulsar Wind Nebulae 
2.1.1. The Crab Nebula 

The Crab Nebula (Ml, NGC 1952) is the first source detected at TeV en- 
ergies and it is considered traditionally the standard candle of the TeV sky. It 



was first detected in 1989 with the Whipple Cherenkov telescope jSOll followed 



by the HEGRA telescopes nlw which reported and confirmed a bright steady 
source (at the instrument's sensitivity) with a flux at 1 TeV of (2.83±0.04)-10^^^ 
TeV^^s^^cm^^ and a spectral index of 2.62±0.02. With the third generation of 
Cherenkov telescopes, namely MAGIC, H.E.S.S. and VERITAS, the Crab Neb- 
ula has been extensively observed, putting constraints to the low and high ends 
of the energy spectrum, size and variability [19, i3i3,|45l]. But even if the Crab 



Nebula is perhaps the best studied/observed source in the TeV sky, it is still one 
of the most mysterious sources. The multi- wavelength nebula spectrum shows 
an undisputed synchrotron nature of the non-thermal radiation from radio to low 
energy y-rays, indicating the existence of relativistic electrons of energies up to 
10^^ eV. The Compton scattering of the same electrons leads to effective VHE 
emission while the synchrotron component is responsible for the radiation from 
radio to relatively low energy y-rays (E<1 GeV) [38]. Recently, the AGILE and 
Fermi collaborations have reported flaring activity above 1 GeV from the nebula 
|3,|45|]. If these flares are related to synchrotron emission of electrons from the 
wind termination shocks, only the highest energy electrons above 100 TeV (in 
fact much higher than 100 TeV) can be responsible for this emission in the E> 100 
MeV band and only these electrons can provide changes on scales of days for any 
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reasonable assumption about the magnetic field. That implies a corresponding IC 
photon flux above at least lOTeV. 

Observations with CTA will be crucial to finally understand the non-thermal 
processes occurring within the Crab nebula. At high energies, CTA should be able 
to constrain the maximum energy of the accelerated electrons and determine the 
cutoff energy with precision. Searches of variation of the spectrum at the highest 
energies should help us to unveil the origin of the Fermi and AGILE detected 
flares at short time scales. These tails may be also visible at the lowest part of the 
CTA-detected energy spectrum if the energy threshold is low enough. 

We evaluated three proposed configurations, CTA_B_ifae, CTA_CJfae and 
CTA_EJfae0 (B, C, E) using CTAmacrosv6 for two different energy spectra. For 
the lowest energies, we used the Crab Nebula spectrum measured by the MAGIC 
telescope. The photon spectrum between 60 GeV and 9 TeV is well described by a 
curved power-law dN/dE=Io(E/300 GeV)[-^-^ '°gi'^(E/300GeV)] ^j^j^ ^ normal- 
ization at 300GeV of Io=(6.0 ± Q2)-lQ-^^TeN-h-^cm-^ , a=-2.31±0.06 and 
b=-0.26±0.07 TeV. To evaluate the cutoff energy, we used the spectrum measured 
by H.E.S.S. between 440 GeV up to 20 TeV. A fit to a power-law function with an 
exponential cutoff dN/dE=Io(E/l TeV)^'^exp(-E/Ecut) results on a differential flux 
normalization at 1 TeV of 0A5±<d.Ql)-lQ-^'^TQW-h-^cm-^, with r=2.39±0.03 
and a cutoff energy Ecut=14.3±2.1 TeV. These two spectral shape are simulated 
and the reconstructed spectra as detected by CTA with the different configura- 
tions are fitted with a MAGIC- and HESS-like function. The results of the fits are 
listed in Table [U 

Fig. |2]a) shows the reconstructed spectrum for the three configurations using 
the MAGIC spectrum. A zoom of the lowest energies is shown on Fig. [2]b). 
The original spectrum is shown in black. The three configurations have similar 
performance although the energy threshold obtained for configuration B and E is 
lower (0.13 TeV) than for C (0.32 TeV), which is optimized to be more sensitive 
at high energies. Between the two first, the balanced configuration (E) seems to 
reconstruct more accurately the original spectrum at low energy (see Table [T]). 
The simulations show that the spectrum at low energy where the IC peak is mea- 
sured should be derived with high accuracy in a relatively short observation time, 
although, configurations B and E provide a lower energy threshold with respect 



Monte Carlo Simulation from IFAE using configuration B - optimized for low energy events, 
C - optimized for high energy events and E - balanced configuration between the lowest and 
highest energies 
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Figure 2: On the left a) MAGIC-like reconstructed spectrum for the Crab Nebula using configura- 
tions B (blue), E (green) and C (red) in 25 h. On the right b) zoom of the lowest energies for the 
MAGIC-like spectral shape 



Conf. B - expected speclmm (CTA], 25.0011. Ecut=14.3 TeV 
Conf. C - expected spectrum (CTA), 25.00li. Ecut=14.3 TeV 
Conf. E - expected spectrum (CTA), 25.00l>. Ecut=14.3 TeV 




100 



energy E(TeV) energy E(TeV) 

Figure 3: On the left a) HESS-like reconstructed spectrum for the Crab Nebula using configura- 
tions B (blue), E (green) and C (red) on 25 h. On the right b) zoom of the highest energies for the 
HESS-like spectral shape 



Conf. B - expected spectrum (CTA), 25.a0h, Ecut=14.3 TeV 
Conf. C - expected spectrum (CTA), 25.aoh, Ecut=14.3 TeV 
Conf. E - expected spectrum (CTA), 25.00h, Ecut=14.3 TeV 



configuration C, allowing a better determination of the IC peak position. 

A similar analysis is shown in Fig. [31 using the HESS-like spectrum shape. 
A zoom of the cutoff region is displayed on the right b) (with fit values listed in 
Table [I]). It shows clearly that similar to the previous case the three configuration 
reconstruct the high energy part of the spectrum correctly, with slight differences 
between them. Simulations with configuration C allows the extension of the spec- 
trum up to 50 TeV, which might be relevant with lower statistics. 

Concerning the variability observed at low energy (up to a few tens of GeV) 
it seems unlikely (given the minimum energy threshold obtained in 25h) to detect 
directly emission for eventual flares like the ones detected by AGILE and Fermi 
LAT. Nevertheless, if those flaring episodes are due to a change in the maximum 
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Table 1: Fit values resulting from the simulation of a HESS-like and MAGIC-like Crab Nebula 



spectrum 
Configuration 




lo (10^11 TeV^icm-^s-i) 




HESS-like 


Ecut 




r 


B 


14.0±0.7 


3.52±0.03 


2.385±0.004 


C 


14.7±0.4 


3.48±0.02 


2.388±0.004 


E 


14.3±0.6 


3.50±0.03 


2.387±0.004 


MAGIC-like 


b 




a 


B 


0.249±0.005 


0.600±0.002 


2.312±0.004 


C 


0.255±0.005 


0.601 ±0.002 


2.307±0.005 


E 


0.260±0.005 


0.601±0.003 


2.308±0.004 



energy reached of the underneath electron population or a change in its spectrum, 
it should be reflected on the IC emission at the highest energies. To evaluate the 
response of the different configuration to this type of effect, we simulate a softer 
cutoff, which translate in a harder spectral shape like dN/dE=Io(E/l TeV)-'^[exp(- 
E/Ecut)]*^'^ for an observation time of 25h. The reconstructed spectrum (see Fig. 
S]) is fit with the corresponding original spectra. 



Fig. |5] shows the 1 and 2a contours for the best-fit parameters. Configuration 
C and B reconstruct the spectral parameters correctly for the two type of functions 
used while configuration E seems to have a worse response when deriving the 
energy cut. Nevertheless, the three configuration lead to compatible results and a 
hardening of the Crab Nebula spectrum should be detectable with any of the three 
proposed configuration in less than 25h. 

Finally we check the capability of CTA to detect changes in the Crab Neb- 
ula size due, for instance, to cooling effects. In that scenario, a enlargement of 
the size (IC produced by low energy electrons around the compact PWN) could 
be expected at very low energy when comparing with the highest energies. We 
used again the B, E and C configurations and considered the case cooling in 
CTAmacrosv6.C with soft (r=2.7) Gaussian (a=0.1) point-like source and a sec- 
ond slightly smaller (a=0.07) Gaussian-shaped source with harder spectral index 
(r=2.0) for an observation time of 25 h. The resulting photon image is then fitted 
with a Gaussian function convolved with the configuration point spread function 
(PSF) and the fit sigma is recovered for two energy bands, a low energy range 
between 0.05 and 0.1 GeV and a high energy range between 1 and 10 TeV. The 
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Conf. B 


expected spectrum (CTA), 25.D0fi, Ecut=14.3 TeV 
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Figure 4: HESS-like reconstructed spectrum for the two type of spectral cutoff considered (see 
text) for configurations B (blue), E (green) and C (red) for a simulated 25h observation . The solid 
line correspond to an exponential cutoff plus power-law spectral shape, whereas the dashed line 
represent a softer exponential cutoff. 




Figure 5: 1 and 2(7 contours for the fit parameters to the measured Crab Nebula using configuration 
B (blue), E (green) and C (red) for the HESS-like spectrum using an exponential energy cutoff 
(solid lines) and a softer one (dashed lines). 
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Table 2: Fit values resulting from cooling effects simulation for a Crab-like source, for 25 h 
observation time. At low energies (LE) the reconstructed Gaussian sigma should be 0.1° whereas 
at high energies (HE) a shrink should be observed ((7=0.07°). 
Conf Ic. u. 5 c. u. 10 c. u. 50 c. u. 100 c. u. PSF 

HE 



B 
C 
E 



0.07 ±0.04 
0.07 ±0.02 
0.07 ±0.03 



0.07±0.01 
0.072±0.006 
0.073±0.009 



0.073±0.007 
0.072±0.004 
0.072±0.005 



0.072±0.003 
0.072±0.002 
0.072±0.002 



0.072±0.002 0.046 
0.072±0.002 0.058 
0.072±0.002 0.058 



LE 



B 
C 
E 









0.14±0.05 





0.10±0.03 
0.1±0.3 
0.2±0.2 



0.11±0.01 
0.1 ±0.1 
0.14±0.08 



0.104±0.007 0.163 
0.10±0.08 0.304 
0.13±0.05 0.304 



results are listed in Table [21 At high energies, the Gaussian is well reconstructed 
for all configurations and Crab Nebula flux level, configuration C being the one 
that reports smaller errors in the fitted sigma. On the contrary, at low energies, the 
flux has to be increased at least up to 10 Crab unit (c. u.) (for 25 h observation 
time) in case of configuration C and E, and up to 5 c. u. for configuration B (as 
expected since B should be optimized for low energy events). 



2.1.2. PWNe Morphology 

To assess the capability of CTA to resolve the morphology of PWNe, differ- 
ent VHE y-ray shapes for given spectral inputs were simulated. The y-ray flux 
is then convolved with the CTA response and the VHE y-ray excess is recovered 
after subtracting the cosmic -ray background. Fig. [6] shows two particular cases. 
On the left, the excess as it will be observed with CTA for configuration E for a 
2-dimension Gaussian-shape source is displayed whereas the right panel shows a 
second source composed by a Gaussian point-like source surrounded by a shell- 
like excess (composite SNR). The reconstructed sigma of different Gaussian- 
shape excesses do not show large uncertainties for the different configurations 
(or different energy bands). Fig. |7] shows one example (1 c. u., 25 h) for the three 
configurations (in green, blue and red for configurations E, C and B respectively) 
in the 1 to lOTeV energy band. 

To evaluate the response of the different configuration for very large sources, 
we simulated a Gaussian source with C7=1.2°, similar to the extension observed 
by H.E.S.S. in Vela X il2|] . The results show a precise reconstructed Gaussian <7 
with errors of the order of 0.001 for all configurations. 
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Figure 6: Gaussian-shape and composite SNR-shape simulations for a source with a Crab-Uke 
spectrum as detected with CTA Configuration E after 50 hours 
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Figure 7: Gaussian extension reconstructed for different configurations (B in red, C in blue and E 
in green) for a 1 c. u. source in 25h 



Finally, the radial profile of a simulated composite SNR is derived for the three 
configurations and two integral fluxes, 1 and 0.1 c. u. in 50 h. For a composite 
SNR at 1 c. u., a fit to double Gaussian as opposite to a single one fits better the 
data for the three configurations (in the whole CTA energy range from 0.05 to 100 
TeV) (Fig. [Uleft) whereas for 0.1 c. u. in 50 h, the components are only clearly 
distinguished for configuration C (in blue in Fig. [8]right). 

In conclusion, no large differences are observed, being the configuration C the 
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Figure 8: Radial profiles derived for a composite SNR simulation (B in red, C in blue and E in 
green). The dashed lines show the fits to a single and a double Gaussian function for each 50h 
simulation. 



most appropriated for composite SNRs. 

2.1.3. Cooling ejfects 

CTA observations were simulated in order to study energy-dependent mor 
phologies in extended sources. The case of HESS J1825-137 is used as an exam- 



ple of morphology shaped by synchrotron cooling. The reference paper is nll\\ . 
from which we take the surface brightness and spectral index distribution. We 
also assume that those distributions can be extended with the same behavior out 
to larger radii than those probed by the H.E.S.S. observation. 

Configuration B and E were used, for which offset-dependent Monte Carlo 
simulations are available. As the majority of the investigated spectra are rather 
steep, we anyway expect these two to be the most important configurations (con- 
figuration D/C are optimized for high energies). 

In order to simulate the spectral evolution and maps, we simulate similar 



wedges to those used in For each wedge, the photon spectrum is a power-law 
(dN/dE = No{E/Eq)-^) with spectral index T. 

Examples of the simulated index distribution are shown in Fig. |9l It can be 
seen that with the assumed index and surface brightness distribution, E will be 
able to probe out to a 2.5° distance with only 5h of observation time. Within 
50h, the error on the index determination is much reduced, but we will not be 
able to probe out to a much larger extension. Configuration B will probe out to 
larger extension and with a better error estimation on the photon index (see Fig. 
[TOl) . This is due to the better sensitivity at lower energies, necessary for the rather 
steep spectra of the outer shells. 

We can also simulate the capability of CTA to disentangle different spectral 
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indices for the same surface brightness distribution. In 20h observation time it 
is possible to distinguish among scenarios of synchrotron cooling or adiabatic 
cooling (constant spectral index). This is shown in Fig. [TT] where x refers to the 
angular distance to pulsar. 



HESS result 
5 hours, CTA-E 
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1-5 



\— HESS result 
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Figure 9: Simulated spectral index as a function of distance from the center of the object in ma- 
genta point. The assumed intrinsic distribution is derived from the H.E.S.S. data (blue) and it is 
shown as a black soUd line. On the left a) 5h of observations with configuration E. On the right b) 
50h of observations with configuration E. 



Simulated CTA subset, 50 hours, CTA-B 
Simulated CTA subset, 50 hours, CTA-E 




Distance (deg) 



Figure 10: Reconstructed spectral index for a 50h observation signal using configuration B vs 
configuration E. Below 2.5° the two distributions are compatible at the 63% level. 



Maps can be derived for a more generalized case of cooling, with an exponen- 
tially decreasing surface brightness distribution and a steepening spectral index 
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ADIABATIC Simulated CTA subset, 20 hours, CTA-E 
SYNCHROTRON Simulated CTA subset, 20 hours, CTA-E 
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Figure 1 1 : Reconstructed spectral index for 20h observation simulated signal using configuration 
E. The two simulated cases, synchrotron or adiabatic cooling show a clear distinguishable spectral 
index distribution 



with increasing distance form the center of the object. The assumed distributions 
are shown in Fig. [121 Dividing the energy range in five logarithmically spaced 
energy bins, we obtain the profiles shown in Fig. [13] (left). The profiles at lower 
energies tend to be more extended for the case of index steepening. In the case of 
constant index distribution, the extension of the source is not energy-dependent. 
An adiabatic cooling scenario would lead to a less extended emission even for the 
same surface brightness distribution (see Fig. [141) . Both configuration E and B 
allow the reconstruction of the expected results. However, due to a better sensitiv- 
ity, configuration E allows us to probe the energy dependent morphology to higher 
energies. 

2.1.4. Population Studies 

To understand the feasibility of population studies, including the complete- 
ness and uniformity of the sample in the Galaxy, we have selected three PWNe 
detected in VHE gamma-rays, namely SNRG21. 5-0.5, HESS J1356-645 and 
Kes 75, which have relatively low 0.3-30 TeV luminosities (of the order of 10^^ 
erg s^^) compared to the other VHE PWNe. Monte-Carlo simulations of these 
three nebulae with CTA have been performed, according to their spectral and mor- 
phological properties as measured with H.E.S.S., for different array configurations 
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Figure 12: Assumed distributions of index (left) and surface brightness (right) for the general case 
of cooling. 
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CTA-E. SO.OOh observation 
Energy 0.05-0.23 TeV and a=0.476+/-0.002 
Energy 0.2-1.0 TeV and a=0.402+/-0.003 
Energy 1.0-4.8 TeV and c=0.376+y-0.005 
Energy 4.8-21.9 TeV and c=0.370+/-0.010 
Energy 21 .9-100.0 TeV and a=0.361+/-0.033 
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Figure 13; Example of cooling scenario profiles. Top: with spectral index change as in Fig. [T2l 
Bottom: adiabatic cooling. 
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Figure 14: The different fit extensions in an adiabatic or in a cooling for the same surface bright- 
ness distribution are shown. In black: The adiabatic case is simulated for a 0.4° extended source. 
In red: A simulated case in which the index changes as in Fig. [T2]in a 1.4° source. 



(CTA_B, CTA_D and CTaJI, assuming an observation time of 20 h with a mean 
zenith angle of 20°. By varying the PWN distances from 1 to 20 kpc, we have 
calculated the so-called horizon of detectability, defined as the distance at which 
the source has a peak significance of 5 C7 (see |16|]). In order to convert these 
horizon estimates into numbers of detectable PWNe with CTA, we have modeled 
the Galactic PSR/PWN distribution according to the [47] logarithmic spiral arms, 
together with the ll26ll galactocentric distribution, and an arm dispersion as a func- 
tion of the galactocentric radius following the Galactic dust model of 113 ill . Any 
potential displacement from pulsar birth place due to the kick velocity has been 
ignored. 

The resulting source distribution is shown in Figure [151 together with the frac- 
tion of visiblell PWNe as a function of the distance to the Sun. 

A large fraction fpwN(~ 0.4-0.8) of SNRG21.5-, HESS J1356- and Kes75- 
like nebulae should be detectable with configuration D and I. de Jager et al, 
2009 have estimated the lifetime of TeV-emitting leptons in such sources to be 



(resp. D) is optimized for low-energy (resp. high-energy) events, while 1 is the balanced 
configuration which provides a better performance over the whole energy range 

^The visibility is defined here as the fraction of the sky seen by an observatory located at the 
same latitude as the H.E.S.S. site, at zenith angles <45°. 
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Figure 15: Fraction of visible PWNe as a function of the distance to the Sun. For each source and 
each CTA configuration considered here, filled and open symbols give respectively the horizons 
of detectability and re solvability, as defined in the text. 



~40 kyr (for B=3/iG, similar to what has been found in several PWNe such as 
Vela X). This gives a total number of detectable G21.5-0.5-, HESSJ1356-645- 
and Kes75-like PWNe with CTA of NpwN ~ 800 fpwN (Tpiv^/40kyr) (Vpwn/^), 
where VpwN is the Galactic VHE-emitting-PWN rate (in units of number of sources 
per century). This implies that ~ (300-600) PWNe should be detected with CTA 
(for configurations CTA-I and CTA-D). These estimates should be taken with 
care, since we have implicitly assumed that the three considered PWNe are rep- 
resentative of the whole population of VHE nebulae, i.e. we have not accounted 
for any time evolution of the PWN sizes and luminosities, which both depend on 
many parameters. 

2.2. Pulsars 

2.2.1. The Crab Pulsar 

The Crab pulsar is the only pulsar detected so far at VHE energies with Cherenkov 
telescopes. In 2008, MAGIC announced the detection of the Crab pulsation above 



25 GeV making use of an novel trigger system [|22|] . Contrary to what one could 
expect from EGRET AGILE [42] (and later on Fermi fl) observations, 
the characteristic double-peaked structure of the Crab light-curve was still visi- 
ble above 25 GeV. Recently, VERITAS Q and MAGIC Hllll] have reported 
the detection of the Crab pulsar at even higher energies, providing measurements 
of the pulsed spectrum which in both cases is well described by a power-law, 
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ruling out in consequence an exponential cutoff shape of the spectrum proposed 
by Fermi LAT. VERITAS spectrum for the total pulsed emission (P1+P2) from 
100 GeV up to 400 GeV follows a power-law with photon index around 3.8 ± 0.5. 
In turn, MAGIC has been able to measure the spectra of each peak separately in 
a broader energy range, from 50 up to 400 GeV. The spectra of both peaks are 
compatible with power-law functions of photon indices 4.0 ±0.8 and 3.42 ± 0.26 
for PI and P2, respectively. 

In order to estimate the CTA potential for detecting the Crab pulsar, we have 
generated a full simulated spectrum for a 50h observation with CTA configura- 
tions B, E and C using CTAmacrosv6. In the case of pulsar observations, the 
signal is expected only at certain pulsar rotational phases, coinciding with the 
peaks seen in the pulsar pulse profiles. One can then make use of the arrival time 
of the photons to further discriminate between y-ray-like events and background 
(both hadronic and nebular). The Crab light curve at VHE energies is character- 
ized by two narrow peaks, extending no more that a 10% of the rotational period 
[[23 . 21]. Therefore we have assumed a 90% background reduction coming from 
a proper timing analysis of the recorded events. Figure[l6]a) shows the results for 
CTA configuration B of simulating the joint fit of the VERITAS and Fermi-LAT 
data reported in [23j, as well as the MAGIC spectrum for the total pulsed emis- 
sion {P\m+P2m) from pll) . It is clear from the figure that CTA should be able to 
collect sufficient photons to reveal the extent of the Crab pulsed emission in en- 
ergy out to at least 1 TeV In fact, the bare detection of the pulsations would take 
less than one hour. Figure [16] b) shows the expected CTA spectra for both Crab 
peaks using MAGIC measurements. Emission of PI will vanish at energies above 
~500GeV. To quantify the goodness of the reconstruction, the simulated spectral 
points were fitted to a power-law function. The results of the fit and its uncertain- 
ties are listed in Table [3] for different CTA configurations. Configurations B and 
E reconstruct the original emission accurately. Configurations without large-size 
telescopes (LSTs), such as C, produces significantly higher errors in the fitting 
parameters. 



Table 3: Fit parameters with their uncertainties to a power-law function (dN/dE ~ Iq x 
(£/100 GeV)-'^) for configurations B, E and C. 



Configuration 


k X 10- 




-V') 




Spectral Index 






(PI + P2) 


PI 


P2 


(PI +P2) 


PI 


P2 


B 


13.62±0.03 


4.34±0.21 


9. 10 ±0.24 


3.57 ±0.03 


3.89 ±0.07 


3.40 ±0.04 


E 


13.61 ±0.04 


4.51 ±0.26 


9.29 ±0.30 


3.58 ±0.04 


3.84±0.10 


3.40 ±0.04 


C 


14.95±1.10 


4.10±1.54 


9.11 ±0.75 


3.71 ±0.08 


3.72±0.85 


3.44 ±0.07 
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Figure 16: Left a): Simulated spectrum of the Crab Pulsar for 50h with CTA configuration B. The 
input models correspond to the joint fit (blue line) of VERITAS and Fermi-LAT data used in lEsll . 
and the MAGIC power-law fit (grey line) from [21]. Right b): Simulated spectra of each Crab 
pulsar peak for 50h with CTA configuration B, using MAGIC power-law fits given in 12 ill . The 
Crab Nebula spectrum from MAGIC is shown in both figures. 



2.2.2. Prospects for new pulsar detections at VHE energies 

With exponential cutoff energies between 0.7 GeV and 7.7 GeV according 
to the 1st Fermi LAT Catalog of Gamma-ray Pulsars [jsl], the current detection 
prospects with CTA for additional Fermi pulsars appear dim at first. Not even the 
Crab pulsar would be detected at VHE energies under this scenario. To illustrate 
this further, Figure [17] shows the spectral fits (power-law with exponential cutoff) 
of the Fermi pulsars taken from [1], in comparison with the standard CTA sen- 
sitivity curve for configuration B in 50 h. The fits of Vela, the Crab pulsar and 
Geminga are indicated explicitly while the shaded area contains the fits for the 
remaining 43 pulsars. 

The Crab detection by MAGIC and VERITAS discussed above, along with the 
fact that at the highest energies measured by Fermi some pulsars already seem to 
deviate from a simple exponential fit, led us to investigate what would happen if all 
the Fermi pulsars were to present a power-law tail like the one seen in Crab. While 
it is unclear that VHE emission is present in every pulsar, CTA could search for 
and place stringent constraints on the putative pulsed emission above a few tens 
of GeV down to a sensitivity of 0.001 c. u. 

To test such Ansatz, we have calculated how may pulsars would be seen with 
CTA depending on the array configuration. In practical terms, we have extended 
the spectra measured by Fermi above their cutoff energy with a power-law tail that 
assumes the same spectral index as the one found for the Crab by VERITAS, when 
a broken power-law is applied to fit both Fermi LAT and VERITAS detections, i.e 
/3 = 3.52. The final extrapolated spectral shapes for 3 out of 46 pulsars are shown 
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Energy [GeV] 

Figure 17: Pulsar spectra fits with power-law and exponential cutoff obtained by Fermi HI] are 
compared to the CTA sensitivity curve for configuration B in 50 h. The fits for Vela, Crab and 
Geminga are indicated explicitly with continuos lines (red, blue and black, respectively). The 
shaded area contains the fits for the remaining 43 pulsars. An extension of the Fermi pulsars' 
spectra to the CTA energy range by adding a power-law tail is shown with long-dashed lines for 
Vela, Crab and Geminga. The photon index /3 of the power-law tail is taken equal to 3.52 (see text 
for explanation). 

in Figure [TTI). Witli such hypotlietical (except for tlie case of the Crab pulsar) 
additional power-law tails all 46 pulsars were then considered as targets for 50 h 
observations with the CTA configurations: B, C and E. We found that 20 pulsars 
would be detected with the configuration B and E; this number reduces to 12 for 
configuration C. This indicates that configurations B and E are better suited for 
pulsar studies than C. Figure [TS] shows how the detectability with configuration 
B depends on the exponential cutoff energy value (as determined by Fermi LAT) 
and the photon flux density at this energy. In conclusion, it seems that under the 
hypothesis of the existence of the VHE Crab-like energy tails, a large fraction 
(up to ~ 40% for configuration B and E) of the brightest Fermi pulsars may be 
detected with CTA. 

On a second step, we can speculate on power-law tails in pulsar spectra with 
other slopes than /3 = 3.52. To test such cases we used broken power-law spectral 
shapes in the form proposed by VERITAS n23ll . The key parameters in this form 
are: Eq - the break energy (usually, quite close to the exponential cutoff energy), 
a -the slope of the photon flux spectrum in the Fermi LAT range well below Eq, 
/3 - the slope of the photon flux spectrum in the tail, i.e. well above Eq (note: 
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Figure 18: Pulsars detectable in 50 h with CTA configuration B under the assumption of Crab-like 
VHE tails as a function of their cutoff energy and flux. 



we define a and /3 with a sign opposite to the original notation in n23 ). For all 
46 pulsars the values of E and a were used based on [3]. With the limit of 50 h 
observation time on one hand and the condition of the overall spectral shape to be 
convex /3 > a on the other hand, the detectability of these pulsars in function of 
/3 was studied for configurations B, C and E. As expected, configuration B is the 
optimal one for all the possible values of /3 ; the second best is configuration E and 
the worst one is C (see Figure [T9l). 

For /3 =3.52 the results presented in the previous paragraph should be identi- 
cal to those of Figure [T8l This is not the case - the latter results give higher per- 
centage of detectable pulsars than the former ones. The reason for this is that each 
analysis was performed by an independent group. The two groups used somewhat 
different analytical functions to approximate the HE- VHE pulsar spectra. More- 
over, their treatment of background effects, of pulsed fraction as a function of 
energy for each pulsar (timing plays an important role in improving the sensitiv- 
ity) etc., was not identical. However, in both cases the detectability of pulsars with 
putative VHE tails is quite high - no less than ~ 40%. Needless to say, there is 
little hope that all gamma-ray pulsars will cooperate in the way described above. 
However, some theoretical models of young and energetic pulsars as well as old 
millisecond pulsars speak in favor of pulsed spectral components located in the 
VHE domain. CTA will be the only facility in near future capable of solving this 
problem, at least to some extent. 
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Figure 19: The detectability of 46 pulsars from |3] in 50 h under the assumption of power-law 
VHE tails as a function of the slope j3 (see text) for configurations B, C and E. 



2.2.3. Globular Clusters 

The two main candidates for detections of unpulsed VHE emission from glob- 
ular clusters are 47 Tuc and Ter 5 due to their large population of MSPs. 

To assess the detectability of these two clusters, we evaluate only the CTA 
response using configuration E and B. Indeed, due to the expected sharp drop of 
the emission at energies of 10 TeV, we assumed that a configuration optimized for 
the high energies would not bring any improvement. Simulations using prediction 
from [25] and [48] (calculated for a population of 100 MSPs) are shown in Fig. [20] 
and[2TJ From the figure it can be derived that the detection of 47 Tuc is expected 
in all the explored models with observation times up to 50 hours. Since the VHE 
spectrum is linearly dependent on the number of MSPs in the system the model 
relative to 47 Tuc is also scaled by a factor of 0.23 (relative to considering only 
the 23 detected MSPs in the system). For this case, CTA will also be able to detect 
it even when high magnetic field values are considered. For the second candidate, 
Ter 5 the fluxes predicted are lower. However, detection can still be achieved with 
a moderately long observation time. For comparison, we show also the best fit 
spectrum from the 90 h observation of HESS J 1747-248 (positionally coincident 
with Ter 5). Different models were re-scaled and discussed in the detection paper. 
Observations with CTA will provide a detailed study of the spectrum and will 
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allow us to associate it or not with the cluster itself (see Fig. [20] and [21] for the 
simulations). Only the response from configuration E and B was investigated. 
The results for the two configurations considered are comparable. Indeed, the low 
energies for which configuration B is optimized are not crucial for the distinction 
of the models investigated here. 

3. Outlook 

The wide range of phenomenological research comprised on the studies of 
PWNe and Pulsars requires diverse specifications. 

For PWN, in general, the simulations described above tend to favor a bal- 
ance configuration (either D or E) in comparison with the configurations B and 
C, optimized for low and high energies, respectively. In particular, the horizon 
of detectability optimizes for configuration D and/or E as well as the sensitivity 
to cooling effects, given the larger range of energy covered with good sensitivity. 
Moreover, although the detailed morphological and spectral investigations show 
an improvement when using configurations which are optimized for low and high 
energies for these particular cases, an intermedium solution seems the most sen- 
sible compromise for PWNe. 

The prospects for high-quality pulsar studies with CTA should be considered 
with caution; the Fermi LAT results are quite meaningful. A sizable progress in 
studies of individual pulsars will be possible with the configurations optimized 
for very low energy range, especially below 50 GeV. None of the planned config- 
urations will be sensitive enough to deal with pulsed radiation of Fermi LAT-type 
energy cutoffs. However, the discovery by VERITAS and MAGIC of pulsed ra- 
diation from the Crab pulsar well above the LAT-inferred exponential cutoff en- 
ergy came as a surprise and yet another example supporting a statement by Trevor 
Weekes twelve years ago that "VHE astronomy is still observation-driven and that 
theoretical VHE astrophysics still lags as a predictive discipline" [491. 

If other LAT pulsars are similar to the Crab pulsar and contain a power-law 
type turnover then configurations B or E may detect a good fraction of these pul- 
sars. More promising results with B or E configurations are expected in future 
studies of globular clusters. The actual performance of B and E configurations 
will depend strongly on (still unknown) population of MSPs in GCs as well as on 
physical properties of ambient fields in the GCs like the strength of magnetic field 
and the field of ambient soft photons. In principle CTA may be able to study both 
spectral shapes as well as angular extension of the VHE radiation from GCs. 
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Figure 20: The expected CTA response in configuration B for the models of 47 Tuc by il25h (green 
hue) and 11481] (black line for high magnetic field in the cluster, red line for low magnetic field). 
The H.E.S.S. upper limit is shown in blue. Top a) the results for assumed 100 MSPs in the cluster. 
Bottom b) rescaled to the number of 23 (as known) MSPs. 
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Figure 21: The expected CTA response for the models of Ter 5 by II25I1 (green line) and ||48 
(black line for high magnetic field in the cluster, red line for low magnetic field). The level of 
HESS J1747-248 is shown in blue. Configuration B is shown on the top a) while the E one is 
shown on the bottom b) 
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